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ABSTRACT 

Two examples are given of models of ion chemistry in reducing atmospheres: 

Titan, which is a satellite of Saturn, and Jupiter, the largest of the gas giants. 

In both ionospheres, layers of hydrocarbon and/or C, H, and N -containing ions 
have been predicted to appear, with larger ions dominating at lower altitudes. 

Altitude profiles are presented for individual Cj- and Cj-hydrocarbon ions and 
larger ions that are represented for example, as C Z H* and C r H f N+. The ac- 
curacy of the predictions is, however, limited by the availability of information 
about the chemistry of these ions. In addition to rate coeffidents and product . 4 
channels for ion-molecule reactions, dissociative recombination coeffidents and 
branching ratios are needed for many hydrocarbon and and related ions. 

1. Introduction 

Solar ionizing photons are absorbed in the thermospheres of planets and satellites, 
and the peak ionization rate is usually found at a column density of about 10 17 cm -2 . 
This altitude may be near or above the homopause, where the transition from trans- 
port by large scale mixing to that by molecular diffusion takes place. On Jupiter the 
peak ionization rate occurs almost 200 km above the homopause, where the major 
constituents are H 2 and H. Although the density of CH 4 eventually exceeds that of 
N 2 at high altitudes on Titan, the ionization peak appears just above the homopause, 
which is near 1050 km. N 2 is the major species at^the homopause, and the mixing 
ratio of CH 4 is 2-10%. .Although hydrocarbon and related ions comprise only a small 
fraction of the ions produced on Jupiter and Titan, models show that such ions may 
be abundant near or below the major ion peaks in both ionospheres. 

The major ions produced in any ionosphere will be transformed by ion-neutral re- 
actions if they are formed in the presence of sufficient densities of neutrals with which 
they can react. In reducing environments, such as the atmospheres of the outer plan- 
ets and their satellites, ionization tends to flow from species whose parent neutrals 
have small proton affinities to species whose parent neutrals have large proton affini- 
ties. On the giant planets, the terminal ions will be large protonated hydrocarbons; 
on Titan, the major species will be both hydrocarbon ions and protonated nitriles and 
amines. In most cases, the available information about ion chemistry is insufficient to 
identify the terminal ions in the bottomside ionospheres, but since larger molecules 
tend to have larger proton affinities, the average size of the ion probably increases 




LOG DENSITIES (cm 3 ) 

Pig. 1. Density profiles of neutral species in the Titan ionosphere. 

with decreasing altitude. In order to improve our understanding of these ionospheres, 
rate coefficients and bran chin g ratios for reactions of hydrocarbons and hydrocarbon 
ions, and neutral and ionic species containing C, H and N are needed, including those 
for dissociative recombination reactions. 

2 . Titan 

In addition to N 2 and CH^, the thermosphere of Titan contains significant densities 
of H2, H, and non-methane hydrocarbons. The background model neutral thermo- 
sphere that we have adopted here is shown in Fig. 1 . The density profiles of N 2 and 
CH 4 are taken from Strobel et ai. (1992), and those of C2H2, C2H4, C 2 H«, and C4H2 
are taken from Yung et al. (1984). The mi-ring ratio of methane assumed in the 
model is 3%. The CO mixing ratio is adopted from the ground-based measurements 
of Gurwell and Muhleman (1995). In our model, we compute the densities of 31 neu- 
tral species and 37 ions, the most important of which are N2, N2H + , H3 , HCNH + , 
Ci- and C 2 - hydrocarbon ions, C4H3 , and two pseudo-ions called C r H^ and C r H v N+, 
which represent, respectively, hydrocarbon ions other than those included explicitly, 
and ions containing C, H, and N (formed, for example, from nitriles and amines). 

We have included photoabsorption and photoionization of N 2) CH 4 , C 2 H 2 , C 2 H 4 , 
C 2 H 6 , N, C, H and H 2 , photodissodation of HCN, and photoionization of CH 3 by 
Lyman alpha. The solar spectrum used is the high solar activity F79050N spec- 
trum of Hinteregger (private communication), and the solar zenith angle is 60°. 
Photoelectron-impact ionization was included for N 2 , CH 4 , H, H 2 , N and C. The 
most important production rate profiles are shown in Fig. 2 . The major ion pro- 
duced is Nj from 600 to 1800 km, and CH 4 above that altitude. Other important 




Pig. 2. Ion production rate profiles in the Titan ionosphere (Fox and Yelle, unpublished). 

ions produced include N + , CH^ } C 2 HJ, and € 2 !!^. 

The ion density profiles were computed including 627 chemical reactions, eddy 
and molecular diffusion for neutrals, and ambipolar diffusion of ions. Most of the 
rate coefficients for the ion-molecule reactions were taken from the compilations of 
Anirich and Huntress (1986) and Anicich (1993). For neutral-neutral reactions, rate 
coefficients were taken from the model of Yung et al. (1984) or from models of the 
interstellar medium (e.g., Millar et al.. 1988). The most important dissociative recom- 
bination coefficients and their sources are shown in Table 1; where no measurements 
are available, values of (3 - 3.5) x 10 -7 (300/T e )°- 5 cm 3 s -1 were assumed. Even where 
dissociative recombination coefficients have been measured, in many cases there is 
little information about the product yields, for which in general there are many pos- 
sibilities, especially for polyatomic ions. Because the ion and neutral chemistries are 
coupled, and because the exothermidty and production of fast atoms depend on the 
channel by which the recombination proceeds, measurements of the product yields 
are important. Here we have assumed that the available channels are populated with 
equal probability. 

The predicted densities of the major ions are shown in Fig. 3. The major ion at 
high altitudes is N 2 H+, and from 1300-1800 km, it is C 2 H 5 . At the peak the dominant 
ion is the pseudo-ion C,H£. The density of C 4 H 3 exceeds that of C,H+near its peak 
at about 750 km. Most other models have predicted that HCNH + would be the 
major ion (Ip, 1990; Keller et al., 1992; Keller and Cravens, 1994). Our model shows 
that HCNH + is lost by reaction with C 4 H 2 , HC 3 N, and NH 3 below 1000 km, and 
by dissociative recombination above that altitude. The most important production 
reactions for C,H+ are reactions of Ci- and C 2 - hydrocarbon ions with acetylene, Cj- 
hydrocarbon ions with methane and a few reactions of Ci- or C 2 - hydrocarbon ions 





Fig. 3. Predicted ion density profiles in the Titan ionosphere (Fox and Yelle, unpublished). 


with ethylene and ethane. Although C*H+ clearly does not represent a single ion, 
the computed peak density of HCNH + is almost an order of magnitude less than the 
total electron density. 

Clearly, a major uncertainty in the ionosphere model is the accuracy of the neutral 
model, including the density profiles computed by Yung et al. (1984). The conclusion 
that HCNH + does not dominate the ion density profile, however, depends only on 
the presence (with a large enough density) of some species that has a higher proton 
affinity than HNC or HCN. or of some other chemical pathway to the protonated 
species. With a fairly thermospheric large mixing ratio of methane, and the presence 
of solar ultraviolet photons, that is all but certain. 

3. Jupiter 

On Jupiter, the major constituents of the thermosphere are H 2 and H, but be- 
low the homopause, which is about 440 km above the ammonia cloud tops, a s m al l 
but significant amount of methane, about 0.001 — 0.01%, is present (e.g., Atreya, 
1986). The models of Kim and Fox (1991; 1994) have shown that, although the 
peak ion production is near 600 km, ultraviolet photons in the wings of the molecu- 
lar hydrogen absorption lines penetrate to below the methane homopanse, producing 
significant ionization of Ci- and Cj-hydrocarbon ions, as shown in Fig. 4. These ions 
are converted to larger hydrocarbon ions with very little loss by recombination. A 
hydrocarbon ion layer forms below the H + and H 3 peaks, in the region 300-400 km 
above the ammonia cloud tops. Altitude profiles of the major hydrocarbon ions at 
local noon are shown in Fig. 5. Insufficient information about the chemistry is known 
to track the evolution of larger hydrocarbon molecules, and two pseudo-ions CaH£ 
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Fig. 4. Ion production rate profiles in the Jovian ionosphere (Kim and Fox^ 1994). 



Fig. 5. Predicted hydrocarbon ion density profiles in the Jovian ionosphere. The H+ and profiles 
peak at higher altitudes and are not shown (Kim and Fox, 1994). 





Table 1. Assumed Dissori^tive Recombination Coefficients 


Species 

Rate Coefficient ( cm 3 s 1 ) 

Reference 

CH + 

1 x 10 _7 (300/T e ) U4i “ 

Mitchell and McGowan (1978) 

CHj 

2.5 x 10 -7 (300/T e )°- 5 

Mul et al. (1981) 

CH+ 

3.5 x 10 _7 (300/T e ) as 

Mul el al. (1981) 

CH+ 

3.5 x 10 _7 (300/T e )°- 5 

Mul et al. (1981) 

CH+ 

11 x io- 7 ( 3 oo/r e ) 0 - 5 

Adams et al. (1984; 1991) 

ct 

3 x 10‘ 7 (300/T e )°- 5 

Mul et al. (1981) 

c 2 h + 

2.7 x 10 _7 (300/T e )°- 5 

Mul and McGowan (1980) 

c 2 h+ 

2.7 x 10- 7 (300/T e )°’ 5 

Mul and McGowan (1980) 

c 2 h+ 

4.5 x 10 _7 (300/T e )°’ 5 

Mul and McGowan (1980) 

C,HJ 

3 x 1(T 7 (300/T e ) 0 - 5 

estimate 

C,HJ 

7.4 x 10- 7 (300/r e )°- 5 

Adams and Smith (1988) 

C,HJ 

3 x 10 -7 (300/T e )° -5 

estimate 

c 2 h+ 

(3.0 - 3.5) x 10“ 7 (300/T e ) 0 5 

estimate 

C 3 H+ 

3.5 x 10 _7 (300/T e )°' 5 

estimate 

C«H+ 

3.5 x 10 -7 (300/T e )°- 5 

estimate 

C,H+ 

3.5 x 10 -7 (300/T e )°- 5 

estimate 

C r H y N+ 

5.0 x 10" 7 (300/T e ) os 

estimate 

HCNH + 

3.5 x 10- 7 C300/T e )°- 5 

Adams and Smith (1988) 

CN + 

3 x 10" 7 (300/T e )°- 5 

estimate 

C 2 N + 

3.0 x 10~ 7 (300/T e )°- 5 

estimate 

c 4 h 3 + 

6.2 x 10- 7 (300/T e )°- 5 

Adams and Smith (1988) 

NH + 

0.43 x 10- 7 (300/T e )°- 5 

Mitchell (1990) 

net 

3 x 10- 7 (300/T e )°- 5 

Mitchell (1990) 

NH 3 

3.3 x 10 -7 (300/T e )°- s 1 

Mitchell (1990) 

NH| 

11 x io- 7 (3oo/t;) 0 - 9 ! 

DuBois et al. (1978) 

N 2 H + 

1.7 x io- 7 (30o/r e )°- 5 

Adams et al. (1984: 1991) 


and C 4 H+ were used to represent ions with three and four C atoms, respectively, and 
subsequent ions in the evolutionary sequence. Although CH, and C 2 H 3 dominate on 
the topside of the layer, at lower altitudes the ions are converted to larger hydrocar- 
bon ions. A more detailed model will require more information about ion-molecule 
and dissociative recombination reactions of many hydrocarbon ions. 
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